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TECHNICAL PAPER 


THE 


EFFECT OF ACCELERATION VERSUS DISPLACEMENT METHODS 
ON STEADY-STATE BOUNDARY FORCES 


I. INTRODUCTION 


When a substructure model is reduced by the Craig-Bampton method a number of 
degrees-of-freedom (DOF) are retained as physical DOF's to provide interface to other 
substructures When more DOF’s are retained in this interface than are actually required, the 
modelT is said to be over constrained. The results of this, when using the displacement 
method are typically an inaccurate distribution of boundary forces. This inaccuracy also 
occurs when there are justifiably many interface DOF’s which result in an indeterminate 
interface. When the acceleration method is used, this inaccuracy is overcome. However, ma y 
people do not fully understand this method and the many ways of implementing it, so 1 
implementation is sometimes haphazard. 

In the space shuttle payload community there has been an increase in the use of over- 
constrained payload models. This has been, mainly, to afford easy recovery of relative 
deflection data between the payload and the shuttle. While there has also been an increase in 
the use of the acceleration method for the recovery of payload displacements and forces, the 
displacement method remains the method used for recovering system displacements and 

forces. 

The purpose of this study is to describe the acceleration method and investigate the 
problem of indeterminate or over-constrained interface substructures. 

This report will first look at the acceleration and displacement methods from the 
example of a simple 2-DOF problem. This will define the methods and describe the 
similarities and differences. Generalizations will then be made to larger systems. A sin p 
two-dimensional, two-beam problem will then be explored and conclusions draw . 
Recommendations for areas of study with space shuttle payload systems will be made. 


II. DYNAMIC BASICS 

Before discussing the acceleration method, it is useful to look at a simple 2-DOF 
system as shown in figure 1. 
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Figure 1. Two-DOF system. 





The equation of motion for this system is: 
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(3a, 3b) 

®, is the frequency of the rigid body mode; i.e., both masses move in unison with no 
deflection in the spnng. o>> is the elastic mode frequency of the two masses moving relative to 

no C dencctioninthefpring SPr '" 8 ' A ngld ' body mode wi " cause n0 internal forces since there is 

be a li„(lr r e fo" Wi " be , made ab0Ut the applied force - F <’>- 11 wil ‘ •* burned to 

»luUo" Therefore? assumpt '0" ® typ.cally made for each time step in a numerical 


F(t) = F s t+F c 

where s stands for slope and c for constant. 

The solutions to the now uncoupled equations in equation (2) are then: 
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(5d) 



The detailed steps for this solution are shown in appendix B. The constants, C i, C 2 ,A, and B 
are determined by initial conditions. 

One can see here that q\ is the generalized displacement of the first mode which is the 
rigid body mode. One can also see that q 2 is the generalized displacement of the second or 
elastic mode. This mode relates to the two masses moving relative to each other, and one can 
see that it has two terms. The first term is the particular or steady-state solution which is 
the response due to the steady-state equilibrium between the applied and inertia forces. The 
second term is the dynamic response due to the harmonic transient. This second term would 
damp out if damping were included in the problem. 

In regard to the generalized accelerations, one can see that q\ is the generalized rigid 
body acceleration due to the applied force and that h is the harmonic transient acceleration 
of the two masses relative to one another. 


Now, the physical displacements and accelerations, *1, *2, *i> and x 2 , can be 
recovered by multiplying the modal displacements and accelerations by the eigenvectors. 
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From these, it can be seen that the displacements, like the generalized displacements 
are made up of rigid-body, steady-state, and transient terms. In fact, the rigid body terms for 
*i and X2 are identical, which they should be since the body is moving as a whole. It can also 
be seen that the steady-state acceleration term is the classic force divided by mass and is 
the same for x\ and x 2 since it is a rigid-body acceleration. 
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These then are the basic equations for the response of a single 2-DOF system. When 
the physical displacements computed this way are used in loads recovery, the method is 
called the displacement method. 


III. ACCELERATION METHOD 


As seen in the previous section, the displacement method uses the eigenvectors to 
compute the physical displacements. Another way to compute the physical displacements is 
to use the physical accelerations computed above and to re-solve the system equations for 
the physical displacements. Reiterating the system equation: 
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Solving for the displacements: 
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Now, inverting 
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is impossible, 


since it is free-free and the determinate is zero, 


however, it is realized that the accelerations were computed free-free and that any time the 
applied forces are in equilibrium with the inertia forces. Thus, one should be able to ground 
the system determinately, invert the stiffness matrix, and not induce any additional forces. 
The 2-DOF system can be grounded determinately at mass 1 by adding a small stiffness (k 2 ) 
to the K ii term of stiffness matrix. The displacements calculated in this way will be called yi 
and y 2 to differentiate the two methods. Doing this, the following is obtained: 
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From this it can be seen that displacement yi is equal to 0 which indicates that the 
ground spring that was added has not deflected, so the assumption of no added forces was 
correct. Also, it can be seen that the rigid body term is gone. This is understandable since 
system was grounded, and it does not affect loads calculations since new forces were not 
added and rigid-body displacements caused no internal forces. One can also see that y 2 1 
exactly equal to the difference between x 2 and calculated from the displacement method. 
For loads calculations, it is the relative displacements which are important so y 2 suits t 
purpose adequately. And finally, y 2 still contains the steady-state and transient terms. 

This is the heart of the acceleration method. The displacements are found in terms ot 
the accelerations. 


IV. EFFECTS OF MODAL TRUNCATION 


At first glance, the benefit of the acceleration method is hidden, since the same results 
have been obtained with much more effort. However, the benefit becomes apparent when 
modes are truncated. If, as above, all modes are used, the two methods are ldentlc ^. 
Suppose the elastic mode is truncated from the 2-DOF problem, from equations (5), (6), and 

(7): 
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(10)’ and h (U) W ** and * 2 ^ Substituted int0 the equations for yi and y 2 (equations (9), 


>1=0, 
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t» u can u 6 seen that Xl and Xl have retained only the rigid-body displacement 

rPflT H / ?h Ch ha ^ e been noted t0 be of little use in loads calculations. However, y 2 has 
retamed the s t ea dy-st ate term which can be a major load contributor. So, it is the accuracy of 
the steady-state displacement term which is the benefit of the acceleration method 


V. GENERALIZATION OF RESULTS 


If one looks closely at the results from the 2-DOF problem in the previous section 
some generalizations can be made to multi-DOF problems. First, ignoring rigid-body 
displacement, the modal displacements are: & & J 
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(14) 


where: 


j = mode number 
i = physical force DOF number 
n = number of physical DOF’s. 
The modal accelerations are: 


9/rigid < tyi j F lit) , 


1=1 


9/eiastic — of A sin {(Ojt) - coj B cos {(Ojt) . 

The detailed steps for these solutions are shown in appendix C. 

Now, the physical displacements are: 
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(16) 



where: 


j = mode number 
i = physical force DOF number 
p = physical displacement DOF number 
m = number of elastic modes. 

The physical accelerations are: 

x P = 1 L ( ppj X <t>‘j F ' (f) 

7=1 L;=i 


+ X [— ^ sin 5 C0S ♦ 

7=1 


( 17 ) 


where: 


j - mode number 
i = physical force DOF number 
p = physical displacement DOF number 
m = number of elastic modes 


r = number of rigid body modes. 

Since it has been determined that it is only the steady-state terms which the 
acceleration method helps, these terms can be isolated as follows: 
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So one can see that the steady-state displacement is determined by the summation 
of contributions from each elastic mode. As modes are truncated, this summation will contain 
more error However, the steady-state acceleration term which causes this displacement is 
dependent only on the summation of rigid-body modal contributions. These rigid-body modes 
are rarely truncated. Therefore, no matter how many elastic modes are truncated, steady- 
state displacements computed in terms of this acceleration will contain no error. Of course the 
transient term of the acceleration is still a summation, so the error introduced to the transient 
displacements by model truncation is unaffected by solution in terms of acceleration. 
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VI. IMPACT ON COMPUTATION OF FORCES 


As has been shown, the truncation of mode shapes of the boundary can lead to errors 
in the displacements computed by the displacement method. The error to the displacement 

vector {*} can be called {e}. Then one can find the error in the forces computed from these 
displacements since: 


{**} = [*] W 


( 20 ) 


so. 


terror * W W • 


( 21 ) 


For a simple 3-DOF problem, if the interface is determinate, then the DOF’s are 
relatively uncoupled: 
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(23) 

So, as shown by these additional terms, the more coupled or indeterminate the system 
the more a force computed for 1 DOF is affected by error in the displacements at other 
DOF s. It should also be noted that because of the nature of the stiffness matrix these 
magnification effects on the displacement error can be quite large, thus producing large force 
errors from seemingly insignificant displacement errors. The acceleration method produces 
extremely accurate displacements which overcomes this problem. 


VII. APPLICATION TO A SIMPLE BEAM 


A somewhat more complicated problem can now be defined for investigation. Take a 
simple two-dimensional beam with element nodes having translational and rotational DOF’s 
(tig. 2). Two of these discrete physical models can be coupled together. Thus, as shown in 
tigures 3 and 4, each model can be divided into interior and boundary DOF’s, coupled together 
at either 2 or 4 translational DOF’s and driven by a force. 

The equations of motion tor these beams and system can be written: 
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1 11 21 31 41 51 61 71 81 



1 11 21 31 41 51 61 71 81 

BERM R 


F(t) 

Figure 3. Two beams coupled at 2 DOF’s. 


BERM B 

1 11 21 31 41 51 61 71 81 



1 11 21 31 41 51 61 71 81 
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Figure 4. Two beams coupled at 4 DOF’s. 
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here: 


i = interior DOF 
b = boundary DOF. 

Equation (26) can then be coupled by a eigenvalue problem. 
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These uncoupled equations can then be solved for {q} and {q}. 

The physical displacements can then be recovered by the displacement or acceleration 
method. Equation (27) is the displacement method, and the acceleration method is shown in 
equation (29). The displacements from the acceleration method will again be called Y to 
differentiate them. 



0 
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F 


[M]\ 

XiB 

x b 


1 XiA l_ 


(29) 


Here, an adjustment must be made so that a comparison can be made between the two 
methods. In the acceleration method, because the stiffness matrix is grounded, the 
displacement of the ground DOF is zero and all other displacements are relative to it. In the 
displacement method, assuming rigid body displacements are zeroed, all displacements are 
relative to the center of gravity. Therefore, a correction is made to the displacement method 
displacements. If the displacement method displacement of the acceleration method grounded 
DOF is distributed by a rigid-body transformation and then subtracted from the displacement 
method displacements, a one-to-one comparison can be made between the two methods. 
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Now, as has already been noted, the acceleration method is only useful in gaining 
accuracy in the displacements caused by the steady-state accelerations. Therefore,, it is 
desirable to isolate this displacement to eliminate confusion with other effects. As seen in the 
acceleration method, it is easy to separate the steady-state, rigid-body accelerations, and 

therefore displacements, from the transients by simply partitioning {</}. This separation is not 
possible in the displacement method. However, it has been shown that if both methods are 
truncated the same, then the transient portions are equal. Therefore, the steady-state and 
transient portions of the displacements can be computed by the acceleration method. This 
transient portion can then be subtracted from the displacement method solution to obtain a 
displacement method steady-state solution. The two methods can then be compared without 
confusion by truncation effects on the transient portions. 


VIII. RESULTS FOR THE SIMPLE BEAM PROBLEM 


For the problem of two discrete physical beam models coupled together as described 
in section VII, the effect of truncation in the system eigenvalue problem, equation (27), was 
investigated. The particular aspect investigated was that of interface forces between the 
substructures or beams as computed by equations (24) and (25). The results were compared 
for both the acceleration and displacement methods and for severely truncated modes (all 
rotational modes truncated) and no modal truncation. These results are tabulated in table 1 . 
Keep in mind that only the steady-state response is compared and that the only truncation is 
at the system eigenvalue level. 

As can be seen, and expected, there is excellent comparison between the two 
methods when all the modes are used. The acceleration method also agrees exactly whether 
modes are truncated or not. However, the displacement method with truncated modes is 
greatly in error. In fact, the error in the applied force on beam A is 90 to 95 percent. 


IX. CRAIG-BAMPTON REDUCTION 


An investigation of what happens if beam B is reduced by the Craig-Bampton method 
is shown in the following. The coordinate transformation for this reduction is: 
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Table 1. Effect of acceleration versus displacement methods on steady-state 

boundary forces. 


[ BEAM A 


ERROR 


REFERENCE 


2 coupled DOF's 

DOF 1 
DOF 41 
DOF 81 

4 coupled DOF's 

DOF 1 
DOF 21 
DOF 41 
DOF 61 
DOF 81 

acceleration 
method 
all modes 

displacement 
method 
all modes 

jl — . — i 

acceleration 

method 

truncated 

displacement 

method 

truncated 

acceleration 
method 
all modes 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000001 

0.00000005 

0.00000001 

0.00000019 

0.00000011 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.02689317 

0.47039463 

0.02689319 

0.00602957 
0.028261 42 
0.44786349 
0.02826125 
0.00602981 

-0.16666667 

0.49999997 

-0.16666665 

-0.00722064 

-0.15944602 

0.50000000 

-0.15944602 

-0.00722064 



BEAM B 


ERROR 



REFERENCE 


abs(abs(REFERENCE)-abs) 


VALUE 


acceleration 

displacement 

acceleration 

displacement 

acceleration 


method 

method 

method 

method 

method 


all modes 

all modes 

truncated 

truncated 

all modes 

2 coupled DOF’s 






DOF 1 

0.00000000 

0.00000000 

0.00000000 

0.01147362 

0.16666667 

DOF 41 

0.00000000 

0.00000004 

0.00000000 

0.00540434 

0.00000000 

DOF 81 

0.00000000 

0.00000001 

0.00000000 

0.01147363 

0.16666665 

4 coupled DOF's 






DOF 1 

0.00000000 

0.00000030 

0.00000000 

0.00057095 

0.00722064 

DOF 21 

0.00000000 

0.00000063 

0.00000000 

0.02842032 

0.15944602 

DOF 41 

0.00000000 

0.00000004 

0.00000000 

0.00320943 

0.00000000 

DOF 61 

0.00000000 

0.00000000 

0.00000000 

0.02841970 

0.15944602 

DOF 81 

0.00000000 

0.00000000 

0.00000000 

0.00057068 

0.00722064 


The coupled system is then: 



QiB 

X b 

XiA 




(32) 


(27), 


I QiB 

The displacements ( X b 

\ xa 

(28), and (29). The interface 


I can 
forces 


be determined by the same methods as equations 
for beam B are computed from equation (31). 


In this study, four different Craig-Bampton boundary conditions were considered. The 
first two were an over-constrained case where all translational DOF’s of beam B were 
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retained as boundary DOF’s; one was for no Craig-Bampton modal truncation and another for 
severely ^uncated Craig-Bampton modes. The third case considered was a Craig-Bampton 
reduction from the physical model down to the coupling DOF’s. Craig-Bampton modes were 
truncated. And finally, the fourth case was a two-step reduction first to the translational 
DOF’s with truncation and then to the boundary DOF’s. 


X. RESULTS OF CRAIG-BAMPTON REDUCTION SCHEMES 


For all four Craig-Bampton cases described in section IX, it was, again, the interlace 
forces that were studied. The acceleration method with and without system modal truncation 
and the displacement method using all modes showed excellent agreement. This data is 
presented in appendix D. 


For all four cases, the only method that indicated error was the displacement method 
using truncated system modes. These results are tabulated in table 2. Keep in m i nd th * 1 this 
data is only for the displacement method using truncated system modes, and that only t e 
steady-state response is considered. The following conclusions can be drawn from this data: 


The severity of truncation of the Craig-Bampton modes does not aftect steady-state 
interface force error. It is the truncation of modes describing the motion across the boundary 
(i.e., system modes) that causes the error. 


Over-constrained Craig-Bampton models, when used with the displacement method, 
introduce sizeable error in the steady-state interface forces. 


Boundary-constrained Craig-Bampton models, when used with the displacement 
method, provide accurate steady-state interface forces. It is noted that they are not as 
accurate as those computed by the acceleration method. 

For boundary-constrained Craig-Bampton models, the error for indeterminate 
interfaces is greater than for determinate interfaces. 

A two-step reduction to the boundary provides just as accurate steady-state interface 
forces as a single-step reduction. 

Finally, a discrete physical model is the limit in over constraint. A discrete physical 
model can be considered to be over constrained by retaining all DOF’s. The error due to 
system mode truncation is substantial since more information is truncated than in Craig- 

Bampton models. 
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Table 2. Effect of Craig-Bampton reduction schemes on steady-state boundary 
forces (truncated displacement method). 


BEAM A 
Physical 

ERROR 

abs(abs(REFERENCE)-abs) 

REFERENCE 

VALUE 


Physical Beam A 



Beam B 

Beam B 

Beam B 

Beam B 



over 

over 

1 -step to 

2-step to 

acceleration 


constrained 

constrained 

boundary 

boundary 

method 

2 coupled DOFs 

no truncation 

truncated 

truncated 

truncated 

all modes 






DOF 1 

0.02689333 

0.02689277 

0.02689245 

0.02689188 

-0.16666667 

DOF 41 

0.47039469 

0.47039462 

0.47039451 

0.47039449 

0.49999997 

DOF 81 

0.02689316 

0.02689277 

0.02689206 

0.02689189 

-0.16666665 

4 coupled DOFs 






DOF 1 

0.00602957 

0.00602840 

0.00602595 

0.00602551 

-0.00722064 

DOF 21 

0.02826125 

0.02826261 

0.02825282 

0.02825380 

-0.15944602 

DOF 41 

0.44786349 

0.44786334 

0.44786317 

0.44786309 

0.50000000 

DOF 61 

0.02826125 

0.02826261 

0.02825274 

0.02825380 

-0.15944602 

DOF 81 

0.00602962 

0.00602892 

0.00602571 

0.00602551 

-0.00722064 


BEAMB 

Craia-Bampton 

ERROR 

abs(abs(REFERENCE)-abs) 

REFERENCE 

VALUE 


Physical Beam A 



Beam B 

Beam B 

Beam B 

Beam B 



over 

over 

1 -step to 

2-step to 

acceleration 


constrained 

constrained 

boundary 

boundary 

method 

2 coupled DOFs 

no truncation 

truncated 

truncated 

truncated 

all modes 






DOF 1 

0.01335352 

0.01335329 

0.00000000 

0.00000000 

0.16666667 

DOF 41 

0.00547892 

0.00547904 

n/a 

n/a 

0.00000000 

DOF 81 

0.01335353 

0.01335329 

0.00000001 

0.00000002 

0.16666665 

4 coupled DOFs 






DOF 1 

0.00046798 

0.00046793 

0.00001097 

0.00001097 

0.00722064 

DOF 21 

0.02862931 

0.02863001 

0.00001101 

0.00001103 

0.15944602 

DOF 41 

0.00327701 

0.00327707 

n/a 

n/a 

0.00000000 

DOF 61 

0.02862883 

0.02862996 

0.00001098 

0.00001099 

0.15944602 

DOF 81 

0.00046812 

0.00046814 

0.00001095 

0.00001097 

0.00722064 
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XI. CONCLUSIONS AND RECOMMENDATIONS 


This study has described the acceleration and displacement methods for use in the 
recovery of coupled-system boundary forces. A simple 2-DOF system has been used for 
illustration. The effect of the choice of method for use with indeterminate or over-constrat e 
boundaries has been investigated. It has specifically looked at results from a simple two- 
dimensional beam problem using both methods. 

In the space shuttle payload community, there has been an increase in the use of 
over-constrained payload models. Much work has been done on the effects of Craig-Bampton 
modal truncation on system displacements and forces, however hide t work hi as I ^ 
system-level modal truncation (i.e., modes across the boundtmy). ^ h " f ' n ^'" 8 h | |1 > v 

indicate that the effect of this system-level truncation is significant. This may be particularly 
true for the 35-Hz system cutoff frequency that is required by the space shuttle. From 1 
study’s findings, the following recommendations can be made: 

The effect of using displacement method recovery for over-constrained payload models 
with the space shuttle directed 35-Hz system cutoff needs to be studied. 

Until this recommended study is conducted, either the acceleration method should be 
used or a secondary Craig-Bampton reduction to only the attaching interfaces should be done 
for payload-coupled loads analyses. 
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APPENDIX A 


Eigenvalues and Eigenvectors of a 2-DOF System 

The homogeneous equation of motion for a 2-DOF system is: 


m i 

0 

Ixx \ k -k 

P')= 

0 ) 

0 

m 2 

\X2 ) [-k k 

1*2 1 

1 0 I 


Assuming a sinusoidal response: 

x\ -A cos (cot-cx) X 2 — A cos ( cot-ot ) 

x x = - co 2 A cos (cot-cx) x 2 = -co 2 A cos ( cot-cx ) 
Substituting into equation (A-l): 


m\ 0 I A I + k —k I A 1 _ 
0 m 2 . i B I -k k \ B 


k-C0 2 m\ -k 


-k k-co 2 m2 J ' ^ 


A _ 


Equation (A-3) has a nontrivial solution only if the determinate is zero, i.e.. 


(A-l) 


(A-2) 


(A-3) 


k ° )2mi k =(k-co 2 mx)(k-co 2 m 1 )-(-k)(-k) = co 2 (co 2 mxm 1 -k(mx+m 1 )) = 0 

-k k-co 2 m 2 /A ^ 


(A-4) 


The roots of equation (A-4) are: 


o 

II 

3 

(A- 5a) 

k(mx+m 2 ) 
2 m i m 2 

(A-5b) 

To define the eigenvectors, the top half ol equation (A-3) is: 


( k-co 2 m\)A - kB = 0 . 


Defining: 


R (k-co 2 mx) _j a2 m± 
A k k 

(A- 6) 

C, = 1-0 = 1 

(A-7a) 

k(mx+m 2 ) mx _ ~ m i 
” m\tn 2 k m 2 

(A-7b) 
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Therefore, the eigenvectors: 


M = 


1 

m 2 


These eigenvectors are then normalized such that <p T M ^> equals unity: 

1 


M 


Inormalized — 


Vwi+m2 

. 1 

[_ VmpFmi 


V5 


W2 


mi(/ni+m 2 ) 

-,/I m 2 — 

m 2 v mi(mi+m 2 ) 




(A-8) 


(A-9) 
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APPENDIX B 


Solution of Uncoupled Equations of Motion of a 2-DOF System 


The uncoupled equations of motion for a 2-DOF system are. 


= , F( f ) -= =(F s t+F c )—=± 


vmi+m 2 


VmpFm 2 


ij 2 +a 2 q2 - (F,t+F c ) Y mi( ” 2 +W2 ) 


The solutions to these equations will take the form: 

q zz C\t 3 +C 2 t 2 +C 3 t+C 4 +Cs sin cot+C(, cos cot , 

q = 3 Cir 2 + 2 C 2 r+C 3 +tuC 5 cos cot—coC(, sin cot , 
q = 6Cif+2C 2 -£tf 2 C 5 sin cot-co 2 C 6 cos cot . 
In equation (B-l) co = 0, therefore: 


6C x t+2C 2 = (F s t+F c ) 


1 


Vmi+m 2 


c,-* 


1 


6 Vmi+m 2 


C 2 = ^ L 


2 V/ni+m 2 


Therefore: 


<7i = 


f_A+f_A + c c 

2 


X 


Vmi+m 2 


Substituting equations (B-3) into equation (B-2): 


n> 2 Cir 3 + a) 2 C 2 r 2 + (m 2 C3+6Ci)t+(« 2 C 4 +2C 2 ) = (F,/+F c ) Y m ,(m!+m 2 ) 


Ci =0 , 

c 2 = 0 , 


(B-l) 

(B-2) 

(B-3a) 

(B-3b) 

(B-3c) 

(B-4) 

(B-5a) 

(B-5b) 

(B-6) 


(B-7) 

(B-8a) 

(B-8b) 
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0) 2 C 3 t+co 2 C 4 = (i F s t+F c ) J — m 2 ( 

V wi(mi+m2) 


(B-9) 


c 3 =i..yC^r: r , 

w 2 V mi(mi+W2) 

c 4 =£t.jzzmzz. ' 

co 2 V mi(mi+m2) 


therefore: 


<?2 = 


TO 
. w 2 


- + 


(A cos C 02 t+B sin <y 2 t) 


V 


mi(mi+/n 2 ) 


Constants and B in equation (B- 11) and constants C 3 and C 4 
determined from the initial conditions. 


in equation (B 


(B-lOa) 

(B-lOb) 

(B- 11) 
-6) are 
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APPENDIX C 


Solution of Uncoupled Equations of Motion of a General System 

The uncoupled equation of motion for a general system is. 

n 

qj+cojqj = ^ <pij Fi(t) , 

1=1 

where: 

j = mode number 
i = physical DOF number 
n = number of physical DOF’s. 

If linear forces are assumed, then the solution to this equation will take the form. 


q = Cir 3 +C 2 r 2 +C 3 f+C 4 +C 5 sin cot+C 6 cos cot , 

(C-2a) 

q — ^C^\t^~i~2F 2 t ~t~F^~i~coF^ cos cot coC ^ sin cot , 

(C-2b) 

q = 6C\t+2C 2 -(0 2 Cs sin cot-co 2 C 6 cos cot . 

(C-2c) 

For rigid-body modes coj = 0, therefore: 


qj = X F i ~ X <Pij( F sitFd) = 6C\t+2C 2 , 
1=1 1=1 

(C-3) 

C\ = tyijFsi » 

6 i=l 

(C-4a) 

O 2 = ~~ y, <t>ijF c i • 
L 1=1 

(C-4b) 

So, for rigid-body motion: 


^gki = *r X 0<7 F « + V ^ 0<7 f c. + C 3 h-C 4 , 
0 ,=1 z 1=1 

(C-5) 

w n n 

frigid = 1 Xr QijFsi + X <PtjFci- 2h QijFiit) ■ 
1=1 «=1 <=1 

(C-6) 
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For elastic-body modes: 


'4j+tfqj = 'Z<t>ijF i . 

i= 1 

Substituting equations (C-2) into equation (C-7): 


(C-7) 


0 )fCit 3 + (l)f C 2 t 2 + ( 0 )f C3+6C\)t + {(oj C 4 + 2 C 2 ) = X QijiFsit + Fd) , 


1=1 

(C-8) 

Ci =0 , 

(C-9a) 

0 

II 

C 

(C-9b) 

n 

0) 2 C 3 t+Q) 2 C 4 = ^ i (pij(F si t+F ci ) , 
1=1 

(C-10) 

i=i 

(C-lla) 

C 4 = ^ X tyijFci • 

°>j 1=1 

(C-llb) 

So, for elastic motion: 


n 

9;«i»uc = - L rX <Ptj F M + A sin (orf) + 5 cos (m,r) , 
°>f i=i 

(C-12) 

Mastic = -0}f A sin ( (Ojt ) - coj B cos (<a/) . 

(C-13) 

Constants A and B in equations (C-12) and (C-13) and constants C 3 and C 4 
(C-5) are determined from the initial conditions. 

in equation 
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APPENDIX D 


Analysis Results 


EFFECT OF ACCELERATION vs DISPLACEMENT METHODS ON 
STEADY STATE BOUNDARY FORCES 

(Beam B - Craig-Bampton - over constrained - no truncation) 


BEAM A 


ERROR 



REFERENCE 

VALUE 

Phvsical 

ahstabs(REFERENCE)-abs) 



acceleration 
method 
all modes 

displacement 
method 
all modes 

acceleration 

method 

truncated 

displacement 

method 

truncated 

acceleration 
method 
all modes 

2 coupled DOFs 

DOF 1 

0.00000000 

0.00000016 

0.00000000 

0.02689333 

-0.16666667 

DOF 41 

0.00000000 

0.00000001 

0.00000000 

0.47039469 

0.49999997 

DOF 81 

0.00000000 

0.00000002 

0.00000000 

0.02689316 

-0.16666665 

4 coupled DOFs 

DOF 1 

0.00000000 

0.00000042 

0.00000000 

0.00602957 

-0.00722064 

DOF 21 

0.00000000 

0.00000000 

0.00000000 

0.02826125 

-0.15944602 

DOF 41 

0.00000000 

0.00000000 

0.00000000 

0.44786349 

0.50000000 

DOF 61 

0.00000000 

0.00000000 

0.00000000 

0.02826125 

-0.15944602 

DOF 81 

0.00000000 

0.00000000 

0.00000000 

0.00602962 

-0.00722064 


BEAMB 

Craia-Bamoton 

ERROR 

abstabs(REFERENCE)-abs) 


REFERENCE 

VALUE 


acceleration 
method 
all modes 

displacement 
method 
all modes 

acceleration 

method 

truncated 

displacement 

method 

truncated 

acceleration 
method 
all modes 

2 coupled DOFs 

DOF 1 

0.00000000 

0.00000001 

0.00000000 

0.01335352 

0.16666667 

DOF 41 

0.00000000 

0.00000003 

0.00000000 

0.00547892 

0.00000000 

DOF 81 

0.00000000 

0.00000002 

0.00000000 

0.01335353 

0.16666665 

4 coupled DOFs 

DOF 1 

0.00000000 

0.00000030 

0.00000000 

0.00046798 

0.00722064 

DOF 21 

0.00000000 

0.00000049 

0.00000000 

0.02862931 

0.15944602 

DOF 41 

0.00000000 

0.00000003 

0.00000000 

0.00327701 

0.00000000 

DOF 61 

0.00000000 

0.00000000 

0.00000000 

0.02862883 

0.15944602 

DOF 81 

0.00000000 

0.00000000 

0.00000000 

0.00046812 

0.00722064 


23 










EFFECT OF ACCELERATION vs DISPLACEMENT METHODS ON 
STEADY STATE BOUNDARY FORCES 

(Beam B - Craig-Bampton - over constrained - truncated) 




ERROR 

bs(REFERENCE 



REFERENCE 

VALUE 

Physical 

abs(a 

:)-abs) 



acceleration 

displacement 

acceleration 

displacement 

acceleration 


method 

method 

method 

method 

method 

2 coupled DOFs 

all modes 

all modes 

truncated 

truncated 

all modes 






DOF 1 

0.00000000 

0.00000000 

0.00000000 

0.02689278 

-0.16666667 

DOF 41 

0.00000000 

0.00000004 

0.00000000 

0.47039459 

0.49999997 

DOF 81 

0.00000000 

0.00000001 

0.00000000 

0.02689277 

-0.16666665 

4 coupled DOFs 






DOF 1 

0.00000000 

0.00000045 

0.00000000 

0.00602840 

-0.00722064 

DOF 21 

0.00000000 

0.00000000 

0.00000000 

0.02826261 

-0.15944602 

DOF 41 

0.00000000 

0.00000000 

0.00000000 

0.44786334 

0.50000000 

DOF 61 

0.00000000 

0.00000000 

0.00000000 

0.02826261 

-0.15944602 

DOF 81 

0.00000000 

0.00000000 

0.00000000 

0.00602892 

-0.00722064 


BEAMB 


~~ ERROR ~ ' “ 

bs(REFERENCE)-abs) 

REFERENCE 

VALUE 

Craig-Bampton 

abs(a 


acceleration 

displacement 

acceleration 

displacement 

acceleration 


method 

method 

method 

method 

method 


all modes 

all modes 

truncated 

truncated 

all modes 

2 coupled DOFs 






DOF 1 

0.00000000 

0.00000000 

0.00000000 

0.01335330 

0.16666667 

DOF 41 

0.00000000 

0.00000005 

0.00000000 

0.00547904 

0.00000000 

DOF 81 

0.00000000 

0.00000000 

0.00000000 

0.01335329 

0.16666665 

4 coupled DOFs 






DOF 1 

0.00000000 

0.00000014 

0.00000000 

0.00046793 

0.00722064 

DOF 21 

0.00000000 

0.00000005 

0.00000000 

0.02863001 

0.15944602 

DOF 41 

0.00000000 

0.00000003 

0.00000000 

0.00327707 

0.00000000 

DOF 61 

0.00000000 

0.00000000 

0.00000000 

0.02862996 

0.15944602 

DOF 81 

0.00000000 

0.00000000 

0.00000000 

0.00046814 

0.00722064 
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EFFECT OF ACCELERATION vs DISPLACEMENT METHODS ON 
STEADY STATE BOUNDARY FORCES 

(Beam B - Craig-Bampton - 1 step to boundary - truncated) 


BEAM A 

Phvciral 

ERROR 

absfabs(REFERENCE)-abs) 

REFERENCE 

VALUE 


1 — 

acceleration 

displacement 

acceleration 

displacement 

acceleration 


method 

method 

method 

method 

method 


all modes 

all modes 

truncated 

truncated 

all modes 

2 coupled DOFs 

DOF 1 

0.00000000 

0.00000043 

0.00000000 

0.02689245 

-0.16666667 

DOF 41 

0.00000000 

0.00000003 

0.00000000 

0.47039448 

0.49999997 

DOF 81 

0.00000000 

0.00000003 

0.00000000 

0.02689207 

-0.16666665 

4 coupled DOFs 

DOF 1 

0.00000000 

0.00000000 

0.00000000 

0.00602595 

-0.00722064 

DOF 21 

0.00000000 

0.00000006 

0.00000000 

0.02825282 

-0.15944602 

DOF 41 

0.00000000 

0.00000000 

0.00000000 

0.44786317 

0.50000000 

DOF 61 

0.00000000 

0.00000000 

0.00000000 

0.02825274 

-0.15944602 

DOF 81 

0.00000000 

0.00000008 

0.00000000 

0.00602571 

-0.00722064 


BEAM_B 

Craig-Bampton 

ERROR 

abstabs(REFERENCE)-abs) 

REFERENCE 

VALUE 



acceleration 
method 
all modes 

displacement 
method 
all modes 

acceleration 

method 

truncated 

displacement 

method 

truncated 

acceleration 
method 
all modes 

2 coupled DOF's 

DOF 1 

0.00000000 



0.00000000 

0.16666667 

DOF 41 

n/a 

n/a 

n/a 

n/a 

0.00000000 

DOF 81 

0.00000000 



0.00000000 

0.1 6666665 

4 coupled DOF's 

DOF 1 

0.00000000 



0.00001097 

0.00722064 

DOF 21 

0.00000000 



0.00001101 

0.15944602 

DOF 41 

n/a 

n/a 

n/a 

n/a 

0.00000000 

DOF 61 

0.00000000 


■ 

0.00001098 

0.15944602 

DOF 81 

0.00000000 

0.00000003 


0.00001095 

0.00722064 
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EFFECT OF ACCELERATION vs DISPLACEMENT METHODS ON 
STEADY STATE BOUNDARY FORCES 

(Beam B - Craig-Bampton - 2 step to boundary - truncated) 


BEAM A 


ERROR 

bs(REFERENCE)-abs) 

REFERENCE 

VALUE 

Physical 

abs(a 

acceleration 
method 
all modes 

displacement 
method 
all modes 

acceleration 

method 

truncated 

displacement 

method 

truncated 

acceleration 
method 
all modes 

2 coupled DOF's 






DOF 1 

0.00000000 

0.00000000 

0.00000000 

0.02689188 

-0.16666667 

DOF 41 

0.00000000 

0.00000004 

0.00000000 

0.47039445 

0.49999997 

DOF 81 

0.00000000 

0.00000003 

0.00000000 

0.02689191 

-0.16666665 

4 coupled DOF's 






DOF 1 

0.00000000 

0.00000000 

0.00000000 

0.00602551 

-0.00722064 

DOF 21 

0.00000000 

0.00000000 

0.00000000 

0.02825380 

-0.15944602 

DOF 41 

0.00000000 

0.00000000 

0.00000000 

0.44786309 

0.50000000 

DOF 61 

0.00000000 

0.00000000 

0.00000000 

0.02825380 

-0.15944602 

DOF 81 

0.00000000 

0.00000000 

0.00000000 

0.00602551 

-0.007220641 


BEAMB 


ERROR 

bs(REFERENCE)-abs) 

REFERENCE 

VALUE 

Craig-Bampton 

abs(a 


acceleration 

displacement 

acceleration 

displacement 

acceleration 


method 

method 

method 

method 

method 


all modes 

all modes 

truncated 

truncated 

all modes 

2 coupled DOFs 






DOF 1 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.16666667 

DOF 41 

n/a 

n/a 

n/a 

n/a 

0.00000000 

DOF 81 

0.00000000 

0.00000000 

0.00000000 

0.00000000 

0.16666665 

4 coupled DOF's 






DOF 1 

0.00000000 

0.00000000 

0.00000000 

0.00001097 

0.00722064 

DOF 21 

0.00000000 

0.00000006 

0.00000000 

0.00001103 

0.15944602 

DOF 41 

n/a 

n/a 

n/a 

n/a 

0.00000000 

DOF 61 

0.00000000 

0.00000002 

0.00000000 

0.00001099 

0.15944602 

DOF 81 

0.00000000 

0.00000000 

0.00000000 

0.00001097 

0.00722064 
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understand this method and the many ways of implementing it, and so its implementation is sometimes haphazard, 
a , J^ stud y Scribes the acceleration and displacement methods for use in the recovery of coupled system boundary forces 
A simple 2-DOF system has been used for illustration. The effect of the choice of method for use with indeterminate or over- 

h “ be “ inV£Sdgated - 11 has looked at results from a simple two^iimensional beam problem using 

. . 1,16 s P ace shuttle Payload community there has been an increase in the use of over-constrained payload models. This has 

been, mainly, to afford easy recovery of relative deflection data between the payload and the shuttle. While there has also been an 
increase in the use of the acceleration method for the recovery of payload displacements and forces, the displacement method remains 
the method used for recovering system displacements and forces. Much work has been done on the effects of Craig-Bampton modal 
truncation on system displacements and forces, however, little work has been done on system level modal truncation (i.e., modes 
across the boundary). The findings of this study indicate the effect of this system level truncation is significant. This may be 
particularly true for the 35-Hz system cutoff frequency that is required by the space shuttle. From this study’s findings, 

— recommendations for areas of study with space shuttle payload systems are made 
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